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ABSTRACT 

Ice clouds influence the climate system by changing the radiation budget and large-scale circulation. 
Therefore, climate models need to have an accurate representation of ice clouds and their radiative effects. In 
this paper, new broadband parameterizations for ice cloud bulk scattering properties are developed for se- 
verely roughened ice particles. The parameterizations are based on a general habit mixture that includes nine 
habits (droxtals, hollow/solid columns, plates, solid/hollow bullet rosettes, aggregate of solid columns, and 
small/large aggregates of plates). The scattering properties for these individual habits incorporate recent 
advances in light-scattering computations. The influence of ice particle surface roughness on the ice cloud 
radiative effect is determined through simulations with the Fu-Liou and the GCM version of the Rapid 
Radiative Transfer Model (RRTMG) codes and the National Center for Atmospheric Research Community 
Atmosphere Model (CAM, version 5.1). The differences in shortwave (SW) and longwave (LW) radiative 
effect at both the top of the atmosphere and the surface are determined for smooth and severely roughened 
ice particles. While the influence of particle roughening on the single-scattering properties is negligible in the 
LW, the results indicate that ice crystal roughness can change the SW forcing locally by more than 10Wm -2 
over a range of effective diameters. The global-averaged SW cloud radiative effect due to ice particle surface 
roughness is estimated to be roughly 1-2 Wm~ 2 . The CAM results indicate that ice particle roughening can 
result in a large regional SW radiative effect and a small but nonnegligible increase in the global LW cloud 
radiative effect. 


1. Introduction 

Ice clouds significantly influence the radiation budget 
of the climate system (Liou 1986; Lynch et al. 2002; 
Baran 2009; and references therein) and also influence 
the large-scale circulations in the atmosphere (Hartmann 
et al. 1992). The quantification of the radiative effects of 
ice clouds has been the objective of numerous studies 
from both satellite-/ground-based observations and nu- 
merical modeling (e.g., Platt et al. 1998; Comstock et al. 
2002; Yang et al. 2007; Lee et al. 2009; Baran 2012; and 
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references therein). However, owing to various uncer- 
tainties, ice clouds still remain one of the least known 
components in the atmospheric system. The uncertainties 
are caused by many factors, from microphysical proper- 
ties, such as the range of particle size distributions (PSD) 
with their associated crystal shapes (habits), to optical 
properties such as extinction coefficient, phase matrix, 
single-scattering albedo, and asymmetry factor. Modeling 
studies by Slingo and Slingo (1988, 1991) show significant 
sensitivity of climate models to the cloud radiative ef- 
fect. Thus, it is imperative to accurately represent both 
ice cloud microphysical and optical properties in atmo- 
spheric numerical models. 

Many ice cloud optical property parameterization 
schemes have been suggested for use in atmospheric 
models. Early schemes parameterized optical properties 
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as functions of temperature based only on limited ob- 
servations (e.g., Platt and Harshvardhan 1988). By as- 
suming ice particles to have spherical shape that is the 
same as liquid clouds (Slingo 1989), the optical proper- 
ties of ice clouds were derived from the Lorenz-Mie 
theory. It soon became apparent that the Lorenz-Mie 
solution was inadequate for representing the scattering 
and absorbing properties of nonspherical ice particles 
because ice cloud particle shapes are complex. Takano 
and Liou (1989) used smooth-faceted hexagonal col- 
umns as the ice crystal model of choice and developed 
a parameterization scheme based on this habit. Ebert 
and Curry (1992) also employed the hexagonal column 
model and parameterized ice cloud optical properties as 
functions of ice water path and the effective particle 
radius. Fu and Liou (1993) developed a parameteriza- 
tion for use in their own radiative transfer code, which 
has subsequently been updated (Fu 1996; Fu et al. 1998). 
In the Community Earth System Model (CESM), which 
uses the general circulation model version of the ra- 
diation code in the Rapid Radiative Transfer Model 
(RRTM) (known as RRTMG), the ice cloud optical 
property parameterization (Mitchell et al. 1996, 2006) 
is based on a combination of the anomalous diffraction 
theory (ADT) approximation and a database by Yang 
et al. (2000). Hong et al. (2009) developed a new pa- 
rameterization for ice cloud bulk optical properties as 
a function of effective particle size for use in climate 
models based on six ice habits and observed ice particle 
size distributions from several field campaigns, again 
assuming smooth particle surfaces. It is important to 
note that the evolution in parameterization scheme de- 
velopment has been toward the use of more ice particle 
habits and more realistic ice particle size distributions, 
but with limited consideration of surface roughness. 
One of the exceptions is Edwards et al. (2007), whose 
parameterization scheme using an ice aggregate model 
accounts for roughened surfaces and was implemented in 
the Third Hadley Centre Atmosphere Model (HadAM3) 
GCM. 

A focus of recent research has been on uncertainty 
reduction in ice cloud properties as more field campaign 
data become available (Heymsfield et al. 2010, 2013) 
from both hemispheres. Additionally, theoretical ad- 
vances help to simulate ice particle optical properties 
more accurately by including more sophisticated ice 
particle shapes and considering ice surface roughness 
(Yang et al. 2000, 2003; Zhang et al. 2004; Yang et al. 
2008; Zhang et al. 2009). How the findings from these 
observational and theoretical advances influence ice 
cloud radiative effect is the focal point of this study, 
specifically the influence of ice particle surface rough- 
ness on radiative transfer and climate simulations. 


Ice particles are most likely not pristine and may 
include air bubbles, internal fractures, and other in- 
homogeneities. In situ observations indicate that ice 
particles tend to have rough surfaces owing to the 
evaporation, sublimation, or riming processes (Cross 
1969; Ono 1969; Davy and Branton 1970; Ulanowski 
et al. 2006, 2012). It is still unknown how ubiquitous the 
surface roughness feature is among ice cloud particles. 
Previous studies have shown the importance of ice par- 
ticle surface roughness effects and the need for them to 
be accounted for in the computation of the ice cloud 
radiative effect (Macke 1993; Yang and Liou 1998). Fu 
(2007) considers the ice surface roughness effect on the 
asymmetry factor and treats it by neglecting the contri- 
bution from delta transmission, obtaining an approxi- 
mate 20 Wm~ 2 impact on the reflected solar flux. A 
recent study by Yang et al. (2012) shows 10-15 W m 2 
differences in upward and downward solar fluxes owing 
to an assumed reduction in the asymmetry factor in- 
duced by ice particle roughness. The primary short- 
coming of the previous studies is an oversimplified 
representation of ice particle surface roughness. An- 
other shortcoming is the application of single-column 
radiative transfer codes to a limited number of indi- 
vidual cases, thereby limiting the representativeness of 
the results. These earlier studies helped, however, es- 
tablish the need to update the effects of ice crystal 
roughness on optical property parameterizations and 
to revisit the radiative forcing estimates. Both issues are 
pursued comprehensively in this paper. We describe the 
data and methodology used in this study in section 2. 
Section 3 shows the simulated ice roughness effects by 
single-column radiative transfer models (RTMs) and 
an AGCM. A summary and conclusions are given in 
section 4. 


2. Data and methodology 

a. Ice particle single-scattering-property database 

In this study, the single-scattering properties for nine 
ice particle habits from the ultraviolet to the far IR are 
provided in a library recently developed by Yang et al. 
(2013). The library is developed using the most up-to- 
date ice refractive index values (Warren and Brandt 
2008) currently available. The habits are those most 
often found in clouds and include droxtals, plates, solid 
and hollow columns, solid and hollow bullet rosettes, an 
aggregate of solid columns, and a small/large aggregate 
of plates. The calculations implement a new treatment 
of forward scattering (Bi et al. 2009) that considers the 
ray-spreading effect more accurately. With this new 
treatment, the delta-transmission term is no longer 
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necessary. What makes the library optimal for use here 
is the availability of three ice particle surface roughness 
conditions for each habit: completely smooth, moder- 
ately roughened, and severely roughened, following the 
approximate approach discussed in Yang et al. (2008). 
The method of particle surface roughness treatment can 
be found in the appendix. The severely roughened ice 
particles compare best with global Cloud- Aerosol Lidar 
with Orthogonal Polarization (CALIOP) backscatter- 
ing data (Baum et al. 2011) and also with Polarization 
and Anisotropy of Reflectances for Atmospheric Sciences 
Coupled with Observations from a Lidar ( PARASOL ) 
data (C.-Labonnote et al. 2001; Cole et al. 2013). Our 
approach in this paper is to perform calculations for the 
completely smooth and the severely roughened cases 
(referred to as “CS case” and “SR case,” respectively), 
and quantify resulting flux differences as the surface 
roughness effect. 


2013). These microphysical data are used for the cal- 
culation of the bulk scattering properties for the in- 
dividual shortwave (SW) and longwave (LW) bands of 
the RTMs. 


c. Calculation of the ice cloud bulk scattering 
properties 

The method to calculate the ice cloud bulk scattering 
properties follows previous studies (Baum et al. 2005, 
2011). The effective particle diameter (Foot 1988; Yang 
et al. 2000; Mitchell 2002) is defined as 


D 


eft 


M 



1 



3 h= 1 

_• 

^min 

2 M 



X 



h = 1 


D . 

min 


V h (D)f h (D)n(D)dD 


A h (D)f h (D)n(D)dD 


( 1 ) 


b. Microphysical data and ice habit mixture 

A total of 14408 particle size distributions (PSDs) 
derived from 11 field campaigns carried out in different 
locations are used in this study. Heymsfield et al. (2013) 
describe the data in more detail. The microphysical data 
were reprocessed where necessary to mitigate the ice 
particle shattering effect on the PSDs. A general habit 
mixture scheme is adopted that uses all nine habits 
(Baum et al. 2011). The same ice habit mixture scheme is 
adopted in remote sensing studies of ice cloud using 
CALIOP (Baum et al. 2011) and PARASOL (Cole et al. 


where D is particle maximum dimension, D mm and D max 
are the minimum and maximum ice particle dimensions, 
V h and A h are the volume and area of an ice particle with 
habit h, Mis the number of habits,/;, is the habit fraction 
of habit h that depends on the ice particle size, and n is 
the number distribution. 

Band-averaged bulk optical properties of ice clouds 
are calculated by integration over both particle size and 
appropriate wavelength band limits. For this study, the 
mass extinction coefficient, the single-scattering albedo, 
and asymmetry factor over a band are provided by 
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Table 1 . Spectral bands in the RRTMG and Fu-Liou RTMs. 
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where A max and A m ; n are the upper and lower wavelength 
limits of a given SW/LW band, (r sca j, is the scattering 
cross section for habit /?, <x ex t,/i is the extinction cross 
section for habit h, Q ex t,h is the extinction efficiency, and 
5(A) is the solar spectrum (Chance and Kurucz 2010) for 
SW bands. For LW bands, the solar spectrum is replaced 
with the Planck function assuming a cloud temperature 
of 233 K. Note that these assumptions are adopted for 
practical reasons because the number of distinct in- 
atmosphere solar spectra and cloud temperatures is too 
large to consider all possible combinations. We find that 
the band-averaged optical properties are largely in- 
sensitive to these assumptions. Two cases of ice particle 
surface roughness (CS and SR) are considered sepa- 
rately. The influence of meteorological conditions (i.e., 
temperature and humidity) on ice particle surface 
roughness is not included in this study because the 
nature of these physical mechanisms is not sufficiently 
known. 

d. The RTMs and AG CM 

The bulk single-scattering properties are parameter- 
ized subsequently as functions of the effective particle 
diameter and implemented in the Fu-Liou and the 
RRTMG RTM codes and the Community Atmosphere 
Model (CAM5.1) AGCM. The Fu-Liou RTM has 6 


solar bands and 12 1R bands; the RRTMG package 
(Iacono et al. 2008) has 14 solar bands and 16 1R bands. 
The spectral bands in the Fu-Liou and RRTMG RTMs 
are given in Table 1. Both use the correlated-/,' distri- 
bution method to treat absorption by atmospheric gases. 

In the next section, results are presented for our new 
ice cloud models to demonstrate the flux differences 
caused solely by assumptions about particle roughness. 

3. Results and discussions 

a. Broadband RTM results 

The differences in ice cloud optical properties be- 
tween the CS and SR cases are shown in Fig. 1. The CS 
results show that the mass extinction coefficient de- 
creases sharply with an increase in D en from 10 to 60 /xm 
and remains under 0.03 nr g 1 for large ice particles. 
The various bands of both RTMs differ little with re- 
gard to the mass extinction coefficient of SR particles. 
As the central wavelength of the spectral band increases, 
the single-scattering albedo gradually decreases from 1 to 
smaller values, and thus becomes more sensitive to D e s. 
The bands with shorter wavelengths have smaller asym- 
metry factors compared with those with longer wave- 
lengths. Not surprisingly, there exist minimal differences 


Asymmetry factor SSA Mass ext coeff (m 2 g 1 ) 


2798 


JOURNAL OF THE ATMOSPHERIC SCIENCES 


Volume 70 





Fig. 1. (top) Mass extinction coefficient, (middle) single-scattering albedo, and (bottom) asymmetry factor as 
a function of effective diameter for (left) the Fu-Liou model solar bands for severely roughened ice particles, (right) 
Corresponding differences between the severely roughened and the smooth ice models. The six bands shown cor- 
respond to the Fu-Liou RTM SW bands in Table 1. 
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Fig. 2. (top) SW, (middle) LW, and (bottom) total cloud radiative forcing (CRF) at solar zenith angle (SZA) = 60° 
as a function of effective diameter and cloud optical thickness as simulated by the Fu-Liou RTM for SR ice particles 
at (left) the TOA and (right) the SFC. 


between the SR and CS values for the mass extinction 
coefficient and the single-scattering albedo, while the 
asymmetry factor shows a much larger effect from 
roughening. Moreover, the SR ice particles have con- 
sistently lower asymmetry factors than the CS particles 
for the same D c ff. The differences between the SR and 
the CS cases generally increase with D^u with the ex- 
ception of bands 5 and 6. Compared with the near-IR 
bands, the shorter wavelength bands (bands 1-4) show 
much larger differences in asymmetry factor between 
the SR and CS cases. This is consistent with previous 
studies showing that the effect of ice particle surface 
roughness is mostly realized as a reduction in the 
asymmetry factor; that is, the CS particles scatter more 


in the forward direction than the SR particles. For the 
mass extinction coefficient, bands 2 and 6 exhibit the 
largest SR-CS difference when D ca is small. Band 6 
also exhibits the largest SR-CS difference in the single- 
scattering albedo when £> e ff is large. 

A test case is designed to illustrate the sensitivity of 
cloud radiative effect to ice particle surface roughness 
in the two RTMs using a common standard midlatitude 
summer atmospheric profile for clear- and cloudy-sky 
simulations. The cloud radiative effect (CRE) is de- 
fined as the difference in the net (SW, LW, and total = 
SW + LW) flux at the top of the atmosphere (TOA) 
and surface (SFC) between the cloudy sky and the clear 
sky: 
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Fig. 3. Ice particle surface roughness effect on (top) SW, (middle) LW. and (bottom) total CRF at SZA = 60° as a 
function of effective diameter and cloud optical thickness for the Fu-Liou RTM at (left) the TOA and (right) 
the SFC. 


CRF = F all sky - F clear = ^cloudy “ F clear) and ( 7 ) 
^all sky = N X ^cloudy + 0- ~ N) X ^dear ’ (®) 

where ivioudy and Tctear are the net fluxes (downward 
minus upward) for cloudy (overcast) sky and clear sky, 
respectively, and N is cloud fraction, assumed to be 1 
here. For this test case, an ice cloud layer is located 
between 12 and 13 km, the solar zenith angle is set to 60°, 
and the surface is assumed to be black for both the SW 
and LW (i.e., surface albedo and emissivity are 0 and 1, 
respectively, at all wavelengths). The CRE is calculated 
for SR and CS cases with both the Fu-Liou and 
RRTMG RTMs as the ice cloud optical thickness ranges 


from 0.01 to 100 and the D a ^ ranges from 10 to 120 gm. 
Figure 2 shows an example of simulated CRE at the 
TOA and SFC from the Fu-Liou RTM as a function 
of D &t f and the ice cloud optical thickness using the 
SR ice particle model. Similar results were obtained 
for RRTMG but are not shown here. The CRE is 
highly sensitive to the changes in ice cloud optical thick- 
ness, while the dependence on D e[[ is weak. These re- 
sults are in agreement with the study by Flong et al. 
(2009). 

The simulated ice particle surface roughness effect 
on CRE (i.e., the difference between the SR and CS 
CRE) for the Fu-Liou and RRTMG RTMs is shown in 
Figs. 3 and 4, respectively. The Fu-Liou and RRTMG 
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Fig. 4. As in Fig. 3, but for the RRTMG RTM. 
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results share many common features. First, both RTMs 
show substantial and dominant roughness effects in the 
SW while minimal influence in the LW. Second, the 
most pronounced roughness effects occur when the op- 
tical thickness ranges between 3 and 12 for all values of 
D e ff. Beyond this higher value of the optical thickness, 
the roughness effect decreases for a given particle size. 
Third, the patterns of the roughness-induced CRE var- 
iation with /J e n and optical thickness from the two 
RTMs are quite similar. Flowever, the Fu-Liou RTM 
tends to have lower roughness-induced SW CRE by 
around 2Wm~ 2 . Moreover, for RRTMG the highest 
value of the roughness-induced CRE center is skewed 
toward a smaller ice optical thickness than for the Fu- 
Liou RTM. These differences could be due to the dif- 
ferent band structures of the two RTMs. 


b. AGCM modeling results 

We now examine the influence of particle roughness 
on the CRE in the CAM5.1 AGCM. Use of an AGCM 
in this analysis demonstrates the impact of surface 
roughness for a broad range of atmospheric profiles 
that occur through the interaction of the RTM in the 
AGCM with other model components. More impor- 
tantly, the AGCM model has the potential to illustrate 
the dynamical/thermodynamical feedbacks that result 
from changing only the ice cloud bulk scattering prop- 
erties. We replace the default ice cloud scattering prop- 
erties in the version of RRTMG used in CAM5.1 with 
our new CS and SR parameterizations. CAM5.1 is run 
for 10 years at 1.9° X 2.5° horizontal resolution and 
with 31 vertical levels. The model is forced with cyclic 
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Fig. 5. (top) Mass extinction coefficient, (middle) single-scattering albedo, and (bottom) asymmetry factor as a function of ice particle 
effective diameter for selected bands in the RRTMG code for (left) the default CAM5 ice cloud model, (middle) the SR ice particle model, 
and (right) the CS model. 


observed sea surface temperature climatology and green- 
house gases to minimize the influence of changing surface 
and atmospheric conditions and mitigate the influence of 
any particular weather events in a given year. The re- 
sults of cloud radiative effect climatology are derived 
from the 10-yr run. 

The optical properties of the new SR and CS ice 
particle models are compared with those of the default 
CAM5 (Mitchell et al. 1996, 2006) in Fig. 5. We find 
considerable differences between the three cases. The 


differences between the CAM5 default scheme and our 
parameterizations are probably not only due to ice 
particle surface roughness, but other factors as well (e.g., 
ice particle habits and size distributions assumed in gen- 
erating the bulk scattering properties). The asymmetry 
factor shows the largest differences; both the new CS or 
SR models display substantial differences from the de- 
fault model currently adopted in CAM5.1. The single- 
scattering albedo of the near-IR bands tends to decrease 
for larger effective diameters. 
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(a) SWCRE (b) LWCRE 




Fig. 6. CAM5.1-simulated 10-yr mean annual total cloud radiative effect (WnL 2 ) and the roughening effects 
(CRE for the severely roughened case minus that of the smooth case) at the TO A: (a) SW cloud radiative effect 
(SWCRE), (b) LW cloud radiative effect (LWCRE), (c) roughening effect on SWCRE, and (d) roughening effect on 
LWCRE. Hatched areas are significant at the 95% significance level. 

Figure 6 shows the CAM5.1 -simulated 10-yr mean 
annual total cloud radiative effect as well as the surface 
roughness effect (again, the difference between the 
SR and CS cases) on total cloud radiative effect at the 
TO A. Flatched areas in Fig. 6 are significant at the 95% 
confidence level. Total cloud radiative effect is shown 
because the CAM5.1 radiation scheme does not treat 
different cloud types individually, thus making it im- 
possible to single out the effect of ice clouds. Table 2 
summarizes the annual global-averaged SW, LW, and 
net cloud radiative effect and the roughness effect. A 
difference of — 1.46 Wm -2 is derived in the net cloud 
radiative effect owing to the roughening of ice particles. 

This is a sizeable value that approaches the magnitude of 
forcing caused by various scenarios of increased C0 2 
concentrations. 

Results from the 10-yr run also indicate that the spa- 
tial and seasonal variations of surface-roughness-induced 
SW cloud radiative effect are substantial. Zonal aver- 
ages (Fig. 7) show that strong SW surface roughness 
effects are found in the tropics and in both hemisphere 
storm-track regions (Fig. 6). The LW differences are 
limited mostly to the tropics. March has the largest SW 
(— 2.35 Wm~ 2 ) and LW (0.67 Wm 2 ) particle roughness 
effects (Fig. 8). The simulated TOA LW CRE shows 


weaker seasonal variations, with a global-averaged value 
of around 25 W m 2 for all seasons. 

The single-column RTM simulations shown earlier 
suggest that roughness has a minimal effect on LW 
CRE. The CAM5.1 simulations, however, suggest some 
prominent LW effects in the tropics. This may be in- 
dicative of either feedback processes or internal vari- 
ability in the AGCM that warrant further exploration. 
Perhaps the perturbation to the SW radiation by ice 
particle roughening induces changes in regional thermo- 
dynamic structure that subsequently lead to persistent 
local circulation anomalies. In turn, these circulation 
changes could lead to further changes in the ice cloud 
properties and, thus, notable LW CRE changes. This 
suggestion is supported by the fact that surface roughness 
leads to both positive and negative anomalies in SW CRE 
in the model (Fig. 6) despite a prominence of SW CRE 


Table 2. Global annual mean of the SR and CS total cloud radiative 
effect as well as the roughness effect (SR minus CS) (W m~ 2 ). 



SR 

CS 

SR - CS 

SW CRE 

-55.67 

-53.84 

-1.83 

LW CRE 

25.34 

24.97 

0.37 

Net CRE 

-30.33 

-28.87 

-1.46 
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Fig. 7. (top) The zonal averages of TOA total CRF and (bottom) the ice particle surface roughness effects (W m 2 ) 
for (left) SW and (right) LW, as simulated by CAM 5.1. 


reduction in the single-column calculations (Fig. 3). On 
a global scale, the roughness effects are smaller but not 
negligible, which suggests that these feedbacks retain 
some influence on global LW CRE. 

The CAM5.1 AGCM simulation results seem to agree 
in general with the offline RTM simulation results but 
are more variable in both sign and magnitude and change 
with time and geographical locations, indicating po- 
tential feedback effects. The CAM5.1 results also pro- 
vide a way of estimating the effects of ice particle surface 
roughening on CRE on a global scale, where the global 
annual average effects on SW and LW CRE are —1.83 
and 0.37 W m 2 (Table 2), respectively. 


4. Conclusions 

While there is an increasing awareness of the exis- 
tence of ice particle surface roughness, few if any studies 
are available that provide an in-depth analysis of its in- 
fluence on cloud radiative effect. This paper attempts to 
provide such a systematic analysis, with new band model 
parameterizations developed using a new database of 
ice particle single-scattering properties from the ultra- 
violet to the far infrared. The new parameterizations 
employ a set of comprehensive microphysical observa- 
tions from field campaigns conducted around the world, 
and adopt ice cloud particle habits similar to those 



123456789 10 11 12 1 
Months 



123456789 10 11 12 1 
Months 



Months 


Fig. 8. (top) The annual cycle of TOA total CRF and (bottom) the ice particle surface roughness effects (W m 2 ) for 
(left) SW and (right) LW, as simulated by CAM 5.1. 
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implemented in satellite remote sensing studies. The 
new band parameterizations are derived assuming that 
all ice particles are either severely roughened (SR) or 
completely smooth (CS). The CS and SR parameteri- 
zations are integrated into the Fu-Liou and RRTMG 
codes and the National Center for Atmospheric Re- 
search (NCAR) Community Atmosphere Model (CAM, 
version 5.1), which uses a version of the latter algorithm. 

The mass extinction coefficient and single-scattering 
albedo are found to be quite similar between the CS and 
SR ice crystal models. The primary effect of ice particle 
surface roughness is in the reduction of the asymmetry 
factor for rough particles compared with smooth parti- 
cles of the same size — a reduction that increases with 
increasing effective diameter D el{ . Thus, consideration 
of the particle surface roughness effect tends to result in 
decreasing the inferred particle size. Both the un- 
certainties associated with ice particle size distributions 
and surface roughness conditions should be considered. 

The simulated cloud radiative effect (CRE) is highly 
sensitive to the variation of ice cloud optical thickness, 
whereas it is only weakly dependent on D efi . The same is 
true for the effect of roughness on CRE. The effect of 
surface roughness in the shortwave (SW) spectrum 
dwarfs that for the longwave (LW) spectrum. Particle 
roughness has the most influence (more than 10 W m 2 at 
SZA = 60°) for cases with D et[ > 80 /jtm and where op- 
tical thickness ranges between 3 and 12. There are dif- 
ferences between the Fu-Liou and RRTMG RTM 
results at various D otl and optical thicknesses that might 
be due to the different spectral bands used in each RTM 
code. 

The AGCM simulations offer a more complete pic- 
ture of the surface roughness effect on CRE than the 
RTM results because they provide spatial/seasonal vari- 
ations and model feedbacks. For example, while the SW 
cloud radiative effect induced by particle roughness is 
still large and dominant, the LW forcing in the tropics is 
larger than expected from the column RTM simulations. 
The apparent LW radiative effect by particle roughening 
in the AGCM simulation may be the outcome of ther- 
modynamic feedbacks that modify cloud properties. This 
is supported by the fact that both positive and negative 
roughening effects on SW CRE occur in the model; we 
will investigate this more fully in future work. The global- 
averaged response due to the adoption of severe particle 
roughening is estimated to be on the order of 1-2 W m 2 
for the SW spectrum but much smaller in the LW. Since 
the global seasonal cloud radiative effect variation is no 
more than 7 W m 2 , the contribution from the treatment 
of ice particle surface roughness is important. A differ- 
ence of — 1.46 Wm~ 2 in the net cloud radiative effect due 
to the roughening of ice particles is small but is similar in 


magnitude to the forcing of greenhouse gases under 
various plausible scenarios. Additional observational 
and numerical studies (e.g., with other GCMs) are re- 
quired to shed more light about the importance of 
crystal roughness on the global radiation budget. 
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APPENDIX 

Representation of Ice Particle Surface Roughness 

The ice particle surface roughness is considered in the 
improved geometric optics method (IGOM; Yang and 
Liou 1998) used for the calculation of the ice optical 
property database (Yang et al. 2013). The surface slope 
of ice particles is represented by a normal distribution 
following Cox and Munk (1954) who studied the optical 
effect of sea surface roughness: 

1 ( 7 2 + 7 2 \ 

P(Z - Z ^ ) = ^ eXP (-^)> (A1) 

where Z x and Z y are the slope variations at the particle 
surface in the x and y directions. By randomly distorting 
the particle surface slope, the degree of roughness is 
controlled by the cr parameter. In this study, we select 
three values of a with the value of 0.5 denoting the se- 
verely roughened particle case, 0.03 denoting the mod- 
erately roughened particle case, and 0 for the completely 
smooth case. Ice particle surface roughness effect is not 
included for particle size parameters smaller than about 
20 where the single-scattering properties are calculated 
by the discrete dipole approximation method. More dis- 
cussion about particle surface roughness can be found 
in Yang and Liou (1998) and Yang et al. (2008). 
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